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Abstract 



J> ' Models where the dark matter component of the universe interacts with the dark energy field 

have been proposed as a solution to the cosmic coincidence problem, since in the attractor regime 



both dark energy and dark matter scale in the same way. In these models the mass of the cold dark 
matter particles is a function of the dark energy field responsible for the present acceleration of the 



universe, and different scenarios can be parameterized by how the mass of the cold dark matter 

Oh; 

particles evolves with time. In this article we study the impact of a constant coupling 8 between 

h ; 

dark energy and dark matter on the determination of a redshift dependent dark energy equation 



of state W£>e{ z ) and on the dark matter density today from SNIa data. We derive an analytical 



^ ■ expression for the luminosity distance in this case. In particular, we show that the presence of 

such a coupling increases the tension between the cosmic microwave background data from the 
analysis of the shift parameter in models with constant wde and SNIa data for realistic values of 
the present dark matter density fraction. Thus, an independent measurement of the present dark 
matter density can place constraints on models with interacting dark energy. 
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I. INTRODUCTION 



Supernovae of type la (SNIa) chronicle the recent expansion history of the universe. 
The accumulated data encrypts information about the composition of the universe and the 
physical properties of its main components, in particular of the dark energy (DE) [1] that 
drives the accelerated expansion today. 

nn n 

Data from SNIa [2|, |3[, the cosmic microwave background (CMB) 4] and large scale struc- 
ture [5y converged to a concordance ACDM model [6], with a nearly flat universe where a 
cosmo logical constant A is the current dominant energy component, accounting for approx- 
imately 74% of the critical density, the remaining being non-relativistic non-baryonic dark 
matter (DM, 22%) and baryonic matter (4%). 

This simple "vanilla" ACDM model is still compatible with data but is not satisfactory 
mainly because it requires a large amount of fine tuning in order to make the cosmological 
constant energy density dominant at recent epochs. 

DE can also be modelled by a scalar field, the so-called quintessence models, either slowly 
rolling towards the minimum of the potential or already trapped in this minimum [l, 7, 8, J 
In this case, the equation of state of DE may vary with cosmological time. 

In spite of the success of the concordance model, one should keep in mind that there exists 
some tension between SNIa and CMB data when the DE equation of state is allowed to be 
different from a cosmological constant. Best fit models for one set of data alone is usually 
ruled out by the other set at a large confidence limit jlO |. SNIa data typically favors large 
values of non-relativistic dark matter abundance Qdm and a phantom-like DE equation of 
state uj < —1. Of course these conclusions are valid only in standard models of DE and DM, 
that is, models where DE and DM are decoupled. This tension has been ameliorated with 
the new data from the Supernova Legacy Survey (SNLS) |3j, as shown in 111 ]. 

An intriguing possibility is that DM particles could interact with the DE field, resulting 
in a time-dependent mass for the DM particles and a modification in its equation of state. In 
this scenario, sometimes called VAMPs (VAriable-Mass Particles) 121 ] . the mass of the DM 



particles evolves according to some 
a linear function of the fielc 
an exponential function 
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12. 
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I7J, 



unction of the dark energy field such as, for example, 
with a inverse power law dark energy potential or 
191 ] with an exponential dark energy potential. Some 
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of these models have a tracker solution, that is, there is a stable attractor regime where the 
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effective equation of state of DE mimics the effective equation of state of DM [15l . [19]. 

The tracker behavior is interesting because once the attractor is reached, the ratio between 
DM energy density p z?jv/ and DE energy density pbe remains constant afterwards. This 
behavior could solve the "cosmic coincidence problem", that is, why are the DE and DM 
energy densities similar today. However, even in these cases a large amount of fine-tuning 



is required for the energy density scale of the scalar potential [20 ]. 

In this paper we investigate the effect of a DE-DM coupling in deriving bounds on the DE 
equation of state and on Qdm from SNIa data. In particular, we show that analyzing the 
SNIa data with a large positive coupling results in a larger value of Q,dmi thus potentially 
increasing the tension with CMB data. This is confirmed by an analysis of the CMB shift 



parameter 



2l| for interacting models with constant equation of state. 



II. A PHENOMENOLOGICAL MODEL 

Variable-mass particles generically arise in models where the quintessence field is coupled 
10 the non-baryonic dark matter field (coupling to baryonic matter is severely restricted 



22j). Such a coupling represents a particularly simple and relatively general form of modified 
gravity: they appear in fact in scalar-tensor models (in the Einstein frame) and in simple 
versions of higher-order gravity theories in which the action is a function of the Ricci scalar. 
From a lagrangian point of view, these couplings could be of the form g(Q)m,oipijj or h(<&)m^(f) 2 
for a fermionic or bosonic dark matter represented by ip and respectively, where the 
functions g and h of the quintessence field $ can in principle be arbitrary. 

Instead of postulating a definite model by choosing two functions defining a DE self- 



coupling potential and DE-DM coupling 23j, one could alternatively follow an approach that 
is more model-independent and closer to observations by introducing a parameterization for 
the DE equation of state W£,e(o) and for the coupling function 5(a), where a(t) is the scale 
factor of the universe. 

The dynamics of the quintessence field, governed by its potential, induces a time depen- 
dence in the mass of dark matter particles. Therefore one would have m = m($(a)) that 
we will parameterize in terms of a function of the scale factor 5(a) in the following way: 

m(a) = m eli S{a ' )dlna ' (1) 



3 



where mo is the particle mass today. In other words, in addition to the usual parametriza- 
tions of the DE variable equation of state we introduce a parametrization for the DM variable 
mass. Just as for the equation of state wde(o), this parametrization allows to study the 
observational data in a systematic fashion in search of new physical phenomena. The phys- 
ical interpretation of the coupling 5(a) is therefore straightforward: it represents the rate 
of change of the DM particle mass, 5 = d In m/d In a. In this paper we will focus on the 
simplest constant coupling 5. 

This variable mass results in the following equation for the evolution of the DM energy 
density p DM - 

Pdm + 3Hp DM - 8(a) Hp DM = (2) 

where H = a/ a is the Hubble parameter. Conservation of the total stress-energy tensor 
them implies that the dark energy density should obey 

p DE + 3Hp DE (l + w DE ) + 5(a)Hp DM = 0. (3) 
Recently, Majerotto et al. p 



(see also |25j) considered the case of constant Wde and 
assumed a tracking behaviour of the DM and DE densities (pde/pdm oc a^) in their analysis 
of SNIa data. In this class of models the modifications of the abundance of mass-varying 
DM particles were studied in [26( and consequences for the evolution of the universe were 
analyzed in 27]. 

However, in quintessence models the equation of state is generally time-dependent. There- 
fore in this work we will study the impact of a constant DM-DE coupling 5 on the determi- 
nation of a redshift dependent DE equation of state wde(z) and on the best-fit value of the 
dark matter abundance today from SNIa data. 

In the case of a constant interaction, eq. (T2J) can be easily solved: 

PD M (a) = p ( g ) M a- 3+5 , (4) 

where p^ M is the non-baryonic DM energy density today. Substituting this solution in eq. 
(J3]) we obtain a differential equation in the scale factor a: 

dpDE 3 (o) _ 4+(5 , . 

—j— + -pDE(l + w D E) + 5p K D J M a =0. (5) 

Before proceeding to an evolving equation of state, it is instructive to study the particular 
case of a constant wde(o) = w, where the solution to eq. ([5]) is given by: 

n ( n \ _ J") -Sq+w) , $ JO) /-3(l+to) n -3+5\ ( R \ 

PDE{a) - p DE a + - - 3w Pdm \ a ~ a ) l b J 
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where p§ E is the DE energy density today The first term of the solution is the usual 
evolution of DE without the coupling to DM. From this solution it is easy to see that one 
must require a positive value of the coupling 5 > in order to have a consistent positive 
value of pDE for earlier epochs of the universe (it is to be noted however that negative values 
of pde could be allowed if the dark energy is in fact a manifestation of modified gravity, see 
e.g. [28]). This feature remains in the case of varying w and we will consider only positive 
values of 8 throughout the paper. 

Furthermore, if 5 < —3w one has a tracking of DE and DM densities at earlier epochs: 

Pdm(o) 5 + 3w 



PdeW 

resulting in 



(7) 



^DE[a) = -— pr — ► - — • (8) 

PdeW + PdmW ow 

Therefore, requiring Qde < 0.1 in the past fixes 5 < — 0.3w in this simple case of constant 
w. It is also interesting that one can analytically compute the scale factor at r where the 
transition from DM to DE occurs in this simple case: 



_ 5 + 3w p§ E 



S+Sw 

(9) 



For instance, for w — — 1 and 5 = 1.6 we find a tr = 0.665 (corresponding to z tr = 0.5). 

For the remaining of this work, we will study an equation of state with the commonly 
used parameterization: 

w D e{z) =w + WiZ, (10) 

in order to compare with other results in the literature. However, it should be stressed that 
SNIa data is not currently very sensitive to time variations in w^e- in some cases, as in the 
comparison with the CMB data, we will use a constant equation of state by fixing W\ = 0. 
We obtain a closed form solution for eq. (jSJ) as a function of the redshift z: 

Pde(z) = p N DE {z) [1 + G(z, w , w u 5)} (11) 

where 

P^)=P { DU 3wiz a + zf 1+W0 - w ^ (12) 
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is the usual evolution of non-interacting (NI) DE density for this parameterization of the 
equation of state and the correction is given by: 

(o) 

Q(z, w , Wl , 5) = 5 e 3wi {3w 1 ) 3(w "- w ^ +s ^T(-3(w - w x ) - 5, 3w 1} 3 Wl (l + z)) (13) 

Pde 



where T(a, Xq, Xi) is the generalized incomplete gamma function: 

>x 



T(a,x ,x 1 )= [ X1 t^e^dt. (14) 

•>xn 



Notice that the correction vanishes in the case of no interaction, 5 = 0, as it should be. 
We have also found solutions for different parameterizations of the dark energy equation of 
state but will not use them in this paper. 

We use the Hubble-free dimensionless luminosity distance dh defined in terms of the usual 
luminosity distance di as 

d h = ^d L . (15) 

The Hubble-free luminosity distance in our model for the parameterization used is given by 
(we are assuming a flat universe throughout the paper): 

d h (z,n%,w , Wl ,S) =(i + z) [ z dz'\ttf\i + z'f + n { S> M (i + z'f- s + 

Jo 1 

(i - - ngkxi + z'f^-^e^' (i + e( z , w , w h 6))} _1 (i 2 6) 

where fi^w anc l ^6° are the dark matter and the baryonic density fractions today, re- 
spectively. In the next section we will compare d^ (or d^ for a given value of the Hubble 
parameter) obtained from our model to observations in order to study the consequences of 
the coupling between DE and DM. 



III. RESULTS 



We will work with two data sets, the so-called gold data set consisting of 157 SNIa [2j and 



the recent SNLS data set of 71 new SNIa with high redshift 



3]. Since these two data sets 



were obtained with different detection techniques and use different methods to extract the 



distance moduli of t 
can be found in 11 



re supernovae, we will analyze them separately. Other recent analysis 



29 



30 



3l|. 



In both cases we used the maximum likelihood method for estimating the cosmological 
parameters. We find the best-fit values by numerically minimizing the likelihood function 



and the 68% confidence level contour plots were obtained by direct integration of the like- 
lihood function. We fixed = 0.05 throughout our analysis. Since the impact of fl^ on 
dh is very weak, the exact value of £1^ is not important. 

The apparent luminosity m of a supernova in terms of its absolute magnitude M is usually 
written as: 

m = M + 51og 10 4 + 51ogc# ~ 1 + 25. (17) 
The absolute magnitude of each SNIa has a spread around a standard value: 



M = M + A 



(18) 



where Mq is the absolute magnitude of a standard SNIa and A represents the corrections 
arising from fits to the color and light curves of each supernova. Hence one can write the 
apparent magnitude as: 

m corr = m — A = 5 logc4 + M (19) 
where the so-called nuisance parameter M is given by 



M(H , M ) = M + 5 log 10 cH, 1 + 25. 
We will use the values of the distance modulus /i given by 



H = m — A — M = 5 log dh 

in order to derive bounds in the parameters in our model. 
For the analysis of the SNIa gold data set, we follow ref. 



(20) 



(21) 



32 



33 



34J ] and estimate the 



likelihood function already marginalized over the nuisance parameter M, which includes Hq, 
for different values of the coupling 5. 

For the analysis of the 71 new SNIa data obtained by SNLS we used the distance moduli 



values Hi obtained from the best fit values (for h = 0.7) o 
and (3) and Mo, and the corresponding errors reported in 
computed as: 



' the correction parameters (a 
3j. The likelihood function is 



C(w ,wi,D, DM ) = exp 



1 ^ (/ij - 5\ogd L (h = 0.7,wo,wi,Q DM ,Zi) - 25) s 



[22) 



i=l " fJ-i 

where a m is the uncertainty associated with the observational techniques in determining 
the magnitudes, a v is associated with the peculiar velocities (and hence negligible for large 
redshifts) and a int is due to the intrinsic dispersion of the absolute magnitudes. 




FIG. 1: Coimtour plots in the wq — £Idm plane for a constant equation of state with 68% confidence 
level for different values of the DM-DE coupling 8 (5 = (dotted line), 5 = 0.2 (dashed line), and 
5 = 0.6 (solid line)). The best fit values are marked by an "*" (5 = 0), "+" {5 = 0.2) and an "X" 
(5 = 0.6), using the gold SNIa (left figure) and SNLS (right figure) data. 

In the case where there is no interaction between DE and DM, we find that the best fit 
values for ^i^ M in a ACDM model (w\ = and Wq = —1) are £lj) M = 0.26 ±0.04 and fi^ = 
0.19 ± 0.02 for the gold set and SNLS data respectively. The best fit from the SNLS data 
is in remarkable agreement with the recent analysis of the 3-year data from the Wilkinson 
Microwave Anisotropy Probe (WMAP3y) 4] alone, which results in = 0.18 ± 0.04, 

whereas the best fit from the gold data set has a somewhat poorer agreement. 

Many parameterizations for the DE equation of state were tested with SNIa data but 
frequently fixing a particular value of or marginalizing over a flat or gaussian prior 
around = 0.27 



10 



32 



33 



34j. However, the agreement between SNIa and CMB 
gets particularly worse for the gold data set when we allow for W£,e{ z ) and VL^ M to vary 
simultaneously without any priors. For instance, fixing w% = 0, that is, a constant DE 
equation of state, we find the best fit values are wq = —2.4 and VIdm = 0-^4 for the gold set 
and tL>o = —1.0 and = 0.21 for the SNLS data. Allowing for an evolving equation of 
state does not significantly alter the fits. A possible tension between SNIa and CMB data 



8 



which was present in the gold data set when one considers models other than the ACDM 



III- 



model practically disappeared in the SNLS data 

We want to investigate in this article the effects of adding a coupling between DE and DM 
on these fits to SNIa data. The coupling will be characterized by the constant 5. We analyze 
first the case of constant w (i.e. w% = 0) and then generalize to a variable equation-of-state. 



A. Constant equation of state 

We show in Figure [T] the 68% confidence level contour plots in the w o x £l$ M plane for 
different values of the coupling 6, keeping w\ = (constant equation of state). Notice that 
the SNLS data results in a better agreement with CMB measurements. One can see the 
existence of a correlation between wq and and that turning on the interaction results 
in a marked tendency towards increasing the best value for fi^jv/, against the CMB results. 

However, a direct comparison of our result for ^l^ M with the CMB results is not possible, 
since the latter has been obtained in the context of uncoupled models. It is therefore 
important to see whether there are upper limits to which are independent of the 

cosmological model. An upper limit to fl^ which does not depend on the background 
cosmology can be obtained from the galaxy cluster dynamics. However the current data 



yield very weak constraints: ref. 



35j gives = 0.30lg;J? so that even Q% = 0.6 is not 
excluded at more than 95% c.l. and actually the strong degeneracy with as allows for even 
higher values. 

On the other hand, for a constant equation of state, we can use the CMB shift parameter 
in order to study the effects of interaction. The shift parameter encapsulates information 
contained in the detailed CMB power spectrum and, more importantly, its measured values 
is weakly dependent of assumptions made about dark e nerg y. It has been used recently to 



put constraints on models with braneworld cosmologies 361 ] . 
The shift parameter R is defined by: 

P Z <iec dz' 

R - a "l W) (23) 

where Zd ec — 1088 is the redshift of matter-radiation decoupling and E(z) = H(z)/H . We 
will use the value R = 1.70 ± 0.03, obtained from WMAP3y data [21] . Clearly for such 
calculations a radiation component with standard conservation equation was appropriately 
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FIG. 2: Countour plots in the wo — £Idm plane for a constant equation of state with 3a confidence 
level for different values of the DM-DE coupling 5. We compare constraints from the gold data set 
SNIa (solid line) and CMB (dashed line) data. 

included. 

In Figs. |2] and [3] we compare the resulting 3a bounds obtained from R and from SNIa 
data. Figure [2] clearly shows the tension between CMB and the gold data set referred to 
earlier. It is clear that the introduction of the coupling makes the CMB and the gold data 
set more incompatible with each other. The coupling favors smaller values of Qdm from 
CMB at the same time favoring larger values of the same quantity from SNIa data. The 
same qualitative behaviour is seen for the SNLS data, although since there is already a good 
agreement with 5 = 0, the situation is less drastic in this case. It is interesting to notice 
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FIG. 3: Countour plots in the wo — £Idm plane for a constant equation of state with 3a confidence 
level for different values of the DM-DE coupling 5. We compare constraints from the SNLS SNIa 
(solid line) and CMB (dashed line) data. 

that assuming a constant Wde — — 1 an d requiring Qdm > 0.1 seems to rule out a coupling 
5 > 0.5 from the bounds arising from the shift parameter. 

B. Variable equation of state 

Allowing wi 7^ does not change the qualitative features of Figure [TJ We show in 
Figure H] the 68% confidence level contour plots in the wq x plane obtained with a 

marginalization over W\ and as expected the allowed region gets somewhat broader, with 
the best fit of fl^ M shifting to larger values. 
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FIG. 4: Countour plots in the wq — £Idm plane marginalized over values of W\ with 68% confidence 
level for different values of the DM-DE coupling 5 (5 = (dotted line), 5 = 0.2 (dashed line), and 
5 = 0.6 (solid line)). The best fit values are marked by an "*" (5 = 0), "+" {5 = 0.2) and an "X" 
(S = 0.6), using the gold SNIa (left figure) and SNLS (right figure) data. 

In Figure [5] we show the marginalized likelihood function over wq and w i in order to 
obtain the 68% confidence level estimation for Q m . The main effect of the coupling is to 
increase the estimation of ^ i s interesting to observe that for the SNLS dataset the 

WMAP3y value for VLdm is rejected at more than 95% c.l. already for 5 > 0.2: however, as 
we anticipated, such a direct comparison is not correct. 

These results are summarized in Table I, where we show the best fit values of wq and 
CIdm with their corresponding la errors obtained by marginalizing over wi and Qdm (in the 
case of Wq) and Wi and w (in the case of VIdm) for both sets of SNIa data. Although the 
best fit values for VL DM are quite large when S is included, the distribution is non-zero even 
for small values of ^dm- 

It is also interesting to study the consequences of the DM-DE coupling in the determi- 
nation of the parameters characterizing the DE equation of state. In Figure [6] we plot the 
68% confidence level contour plots in the Wq x wi plane obtained by marginalizing the like- 
lihood function over VL DM using a gaussian prior VL DM = 0.18 ± 0.04, as obtained from the 
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TABLE I: Best fits obtained for wo and Odm for both sets of SNIa data. 
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FIG. 5: Plot of the likelihood function for Qdm with wq and u;i marginalized for different values 
of the DM-DE coupling 5 (5 = (dotted line), 5 = 0.2 (dashed line), and 5 = 0.6 (solid line)) using 
the gold SNIa (left figure) and SNLS (right figure) data. 



WMAP3y data. The effect of the coupling is very small when the prior on the DM density 
is taken into account. 



IV. CONCLUSIONS 



There is a vast amount of work studying the possibility of having an interaction between 
the dark energy and the dark matter components of our universe. In this paper we analyzed 
a simple model for the interaction between these two fluids, in which the mass of the dark 
matter particles increases at a constant rate. We have shown that introducing a coupling 
between the dark energy component of the universe with dark matter particles has the effect 
of increasing the best fit values of the DM density today obtained from current SNIa data. 

We performed a comparison with CMB data for the case of constant equation of state, 
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FIG. 6: Countour plots in the wo — w\ plane marginalized with a gaussian prior over values of 
Qdm with 68% confidence level for different values of the DM-DE coupling 5 (5 = (dotted line), 
5 = 0.2 (dashed line), and 5 = 0.6 (solid line)). The best fit values are marked by an "*" (5 = 0), 
"+" (<5 = 0.2) and an "X" (5 = 0.6), using the gold SNIa (left figure) and SNLS (right figure) data. 

using the shift parameter. We found that the introduction of the coupling results in a poorer 
compatibility between CMB and SNIa data. Our results showed that assuming a constant 
wjje — —1 and requiring VL DM > 0.1, a coupling 5 > 0.5 seems to be ruled out. This must 
be checked by model-independent measurements of Om from large scale structure and also 
with a more careful analysis of the CMB data including more observables in addition to the 
shift parameter. We also found that introducing the coupling does not change significantly 
the determination of the DE equation of state when a prior on Qdm is adopted. 

We worked with a simple parameterization of the dark energy equation of state and 
for the DM-DE coupling but we believe that our results are fairly general for the type of 
interaction that we introduced. It would be interesting to extend our analysis to more 
general parameterizations for both the interaction and the equation of state available in the 
literature. However, it should be stressed that SNIa data is not currently very sensitive to 
a time- varying wde- 
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